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A model of postphotodissociative monomolecular (gemi- 
nate) recombination of heme proteins with small ligands (NO, 
O2 or CO) is represented. The non-exponential decay with 
time for the probability to find a heme in unbound state is 
interpreted in terms of diffusion-like migration of ligand be- 
tween protein cavities. The temporal behavior for the proba- 
bility is obtained from numerical simulation and specified by 
two parameters: the time r re (, of heme-ligand rebinding for 
the ligand localized inside the heme pocket and the time r esc 
of ligand escape from the pocket. The model is applied in 
the analysis of available experimental data for geminate re- 
oxygenation of human hemoglobin HbA. Our simulation is in 
good agreement with the measurements. The analysis shows 
that the variation in pH of the solution (6.0 < pH < 9.4) 
results in considerable changes for r re b from 0.36 ns (at 
pH = 8.5) up to 0.5 ns (pH — 6.0) but effects slightly on 
the time r esc (r esc w 0.88 ns). 



I. INTRODUCTION 

The binding reactions between myoglobin (Mb) or 
hemoglobin (Hb) and small ligands (NO, O2 or CO) are 
objects of extensive investigations during many decades 
because of the great functional importance of the heme 
proteins for living systems [1-3]. In the investigations 
a special attention is paid to the heme-ligand recombi- 
nation process going after fast photodissociative bond 
breaking between a ligand molecule and an ion Fe + + 
located in the center of a heme (Fe-protoporphyrin IX 
complex). The kinetic study of the postphotodissocia- 
tive recombination allows to obtain detailed information 
on the protein-ligand interaction mechanism, the protein 
structure, the allosteric effect and the medium influence 
on the recombination efficiency (see, for example, refer- 
ences [4-9]). 

The heme is well- wrapped in protein helixes, which 
prevent the iron from the solvent and hinder the ligand 
migration through protein matrix. On a sufficiently long 
time scale (at t < 100 ns in the case of Hb) after dissocia- 
tion, when the ligand is not managed to leave the protein 
and to move significantly away from the parent heme, the 
recombination is a monomolecular reaction designated 
usually as a geminate recombination (GR) [5] . Schemat- 
ically, the GR can be written as [10]: 
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where A is the bound heme-ligand state. The substates 
B and {Ci, . . . , C n } form the unbound state. Each of 
these substates corresponds to ligand localization in an 
individual cavity of protein. The substate B answers to 
the residence of ligand inside the heme pocket (the cavity 
nearest to the iron on the distal side of heme) . The rate 
constants k re b an d k esc specify two competing processes: 
the irreversible heme-ligand rebinding for the ligand lo- 
calized inside the heme pocket (that is, the transition 
from the substate B to the state A) and the migration of 
unbound ligand between the heme pocket and other pro- 
tein cavities (the transitions between B and the substates 
{Ci, ... ,C n }). Immediately after photodissociation the 
unbound ligand is in the substate B. Therefore the quan- 
tity k esc can be associated with ligand escape from the 
pocket. In general, the GR is essentially determined by 
the specificities of heme-ligand interaction (including the 
spin restriction effect, the position and the orientation 
of ligand with respect to the heme plane) [5,11-13] and 
the effect of residues surrounding the heme [14-19]. Im- 
portant factors for the heme-ligand rebinding are also 
the state of tertiary [2,20-22] and quaternary [2,23-25] 
structures of protein, the conformation transitions in pro- 
tein [2,10,26-28] and the solvent impact [23,29-31]. As 
a consequence, the kinetic curve (that is, the probabil- 
ity P (t) to find the heme in unbound state) of GR is 
a non-exponcntially decaying function of time [32-35]. 
After realization of the geminate stage a portion P s of 
the hemes remains in unbound state: P s < 0.01 for NO, 
P s ~ 0.1 0.2 for 2 and P s ~ 0.5 1.0 for CO. The 
quantity P s characterizes the efficiency of ligand escape 
from the protein to the solvent. 

Molecular dynamics simulations [11,36-40] show that 
the movement of unbound NO, O2 or CO ligands in 
heme protein can be associated both with ligand trap- 
ping for a significant time in individual cavities and with 
rare jump-like transitions between adjacent cavities. It 
implies a fast establishment of equilibrium for the prob- 
ability distribution of ligand within individual cavities. 
The establishment occurs on a time scale comparable to 
the mean time interval t w between the collisions of ligand 
with cavity walls. At room temperatures the time t w lies 
in the subpicosecond range (t w ~ 0.1 ps for NO in the 
heme pocket of Mb [41]). The ligand redistribution be- 
tween protein cavities is observed on a longer time scale 
ranging from several tens of picoseconds (~ 40 ps for NO 
in Mb [41]) up to several tens of nanoseconds (~ 50 ns for 
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CO in Hb [42]). Unfortunately, in practice the detailed 
molecular dynamics simulation can not be implemented 
to the GR due to enormous computational efforts. 

In the study we apply an alternative approach based on 
the diffusion approximation to ligand migration in pro- 
tein. Such an approximation is valid for times t ^> t w 
when the deterministic nature of ligand motion can be 
ignored. Here the interval t w can be recognized as a 
correlation time. The diffusion-like character of ligand 
migration in the heme proteins can be a reason of the 
non-exponential temporal dependence for the probability 
P(t) [43-45]. For instance, a two-dimensional diffusion 
is demonstrated for CO in Mb [45] to explain the power- 
law kinetics to be observed in the experiment. Generally, 
reaction (1) can be represented in a three-dimensional 
diffusion approximation by equation 

(JTi 

— = V (£>Vn) - R reh n (2) 

with the diffusion coefficient D = D(x,y,z). The quan- 
tity n — n(x,y,z,t) is the probability density of un- 
bound ligand in the protein. The stepwise function 
R re b — Rreb{x,y, z) specifies the heme-ligand rebinding 
and equals to k re b inside the heme pocket or to zero oth- 
erwise. 

In order to solve diffusion equation (2) and to follow 
the evolution of GR we use a simple model proposed re- 
cently in [46] . The model reproduces dynamics of random 
walk of particle in porous media (such, for instance, as 
glass- like matrices [47-49] ) and takes into account an ini- 
tial retention of ligand inside the heme pocket (that is, in 
the substate B). In the absence of heme-ligand rebinding 
the substate B is realized at times t < r esc (r esc = 1/ k esc 
is the time of ligand escape from the heme pocket to 
others cavities). Only on a longer time scale (t > r esc ) 
the ligand succeeds to leave the pocket and to migrate 
over protein cavities. Due to the diffusion nature of the 
migration the time T esc can be specified in terms of the 
diffusion coefficient D. 

The approach is implemented with the help of a numer- 
ical simulation where the unbound ligand is represented 
by a structureless particle. For simplicity, in the simula- 
tion we make some assumptions. The ligand migration is 
assumed to be restricted to the distal side of heme. The 
ligand motion (realized on a short time scale t < t w ) 
inside the heme pocket is represented by a unforced dis- 
placement of the particle within a restricted hemispheric 
region of space. At t„ < t < r esc the ligand trajec- 
tories are effectively mixed in the configurational space, 
resulting in a homogeneous distribution for the ligand 
inside the cavities. Hence, the probability of irreversible 
heme-ligand rebinding is accepted to be uniform for the 
whole heme pocket. We take into account also that on 
the time scale t 3> t w the fast intracavity displacements 
of ligand for the substates {C\, . . . , C n } do not influence 
essentially on the GR kinetics and can be ignored in the 



simulation. Therefore the ligand displacement exterior 
to the heme pocket is simulated as a random walk (that 
is, as a Brownian-like motion) of the particle outside the 
hemispheric region. This walk is a spatially homogeneous 
diffusion with the diffusion coefficient D. We neglect also 
the structural transformations (such as a shift of the iron 
with respect to the porphyrin ring plane) at the confor- 
mational transition of protein between the unliganded 
and liganded states. According to the model, the tem- 
poral behavior for the probability P(t) can be specified 
in terms of two parameters: the time r esc and the time 
T re b = l/k r eb of heme-ligand rebinding. The description 
of the model is represented in Section 2. 

In order to demonstrate the usefulness of such an ap- 
proach to the GR of heme proteins we apply the model to 
the analysis of available experimental data. We analyze 
the measured recombination kinetics and the efficiency 
for a postphotodissociative GR of human hemoglobin 
HbA [50,51]. These measurements were carried out at 
various pH values of the solution. Here we determine the 
times T re b and T esc as functions of pH and estimate the 
influence of solution properties on the heme-oxygen re- 
binding, the migration of oxygen molecule in hemoglobin 
and the efficiency of oxygen escape from the protein. The 
association of the times r re b and r esc with the time of a 
bimolecular recombination process for hemoglobin is an- 
alyzed. The results of simulation and their analysis are 
represented in Section 3. 

II. DIFFUSION-BASED MODEL OF GEMINATE 
RECOMBINATION 

The movement of unbound ligand is considered in a 
Cartesian coordinate system xyz attached rigidly to the 
heme group of atoms. The system origin is superposed 
on an iron atom located in the middle of heme porphyrin 
ring. The x and y axes are aligned with the heme plane. 
The positive direction for the z axis corresponds to the 
distal side of heme. The ligand migration in protein is 
simulated as a probability redistribution for the ensemble 
of structureless particles over a three-dimensional hemi- 
spheric space with z > 0. As in [52], in our simulation 
the heme pocket is represented by a hemispheric region 
(designated here ge) of radius p. At an initial 

time instant the particle is uniformly distributed inside 
the cage. 

The individual particle to be exposed to a sequence of 
^-shaped uncorrelated kicks executes a random walk in 
the space. As for the Brownian particle, each kick results 
in an abrupt change in the particle velocity. Between the 
kicks the particle is in unforced motion. On a time inter- 
val Atk = tk+i — tk (tk is the time instant of action for 
fc-th kick) between adjacent kicks the particle is specified 
by the velocity and the length Lk of free path (note 
that Atk = Lk/ \vk\ )• Then the radius vector r(tk+i) of 
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particle for the time point of fc+l-th kick can be obtained 
from iteration procedure 



r (tfc+i) = r(t fe ) + 



Kl 



(3) 



where the radius vector r(tfc)is given for the time instant 
of fc-th kick. The projections Vj^ {j — x, y, z) of the 
velocity onto the coordinate axes and the length Lk 
are accepted to be independent random quantities, new 
values of which are generated at each kick. The quantities 
Vj.k is obtained from the Maxwell distribution 



Pm (vj,k) 
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mv 



2nkT 



exp 
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(4) 



Here to is the particle mass and T is a protein tempera- 
ture. At an attainment of the z = plane bounding the 
space, a new particle velocity with v z > is regenerated 
in accordance with distribution (4) . 

The choice of free path length is dictated by the parti- 
cle location in the space. Within the hemispheric cage the 
particle displacement is unforced and the particle under- 
goes no kicks. The length Lk is determined then from the 
ballistic trajectory of particle between the cage bound- 
aries. In this case the length is comparable to the cage 
size p. We accept here that the mean time Th — (Ath,k), 
during which the particle crosses the cage, can be asso- 
ciated with the time interval t w between the collisions of 
ligand with heme pocket walls: Th ~ t w . 

Exterior to the cage, the particle is exposed to uncorre- 
lated kicks. The absence of correlation between the kicks 
implies that the quantity Lk is distributed according to 
the exponential law: 



P{L k ) = ^cxp ( — -j- 



(5) 



where A = (Lk) is the mean length of free path for the 
particle displacement outside of the cage. The mean time 
t c between adjacent kicks and the length A are related to 
the diffusion coefficient D = (LD/6t c by equations: 



T c = 



6mD 
irkT 



X = 3D 



2m 
~vkT 



(6) 



(7) 



Thus, the spatial displacement of the particle is ob- 
tained from iterative equation of motion (3) and depends 
on the random sampling of variables v x ±, v y ^, v Zj k and 
L fe , the statistical distributions for which are specified by 
the parameters m/T, D and p. As mentioned above, un- 
der the conditions typical for the heme proteins (that is, 
the temperature, the ligand mass and the distinctive sizes 
of heme pocket) the times r c and th to be accepted here 



as correlation times are negligibly short as compared to 
the characteristic times of GR. The length A is essentially 
small as against the size p of hemispheric cage. Hence, 
the temporal behavior for the probability redistribution 
of ligand in heme protein can be described in terms of 
the diffusion-based approach. 

Our model reproduces dynamics of ligand migration 
over protein cavities. Initially, the ligand is retained in- 
side the heme pocket and the root-mean-square displace- 
ment S(t) = \J (|r(t) — r(0)| 2 ) of ligand from the initial 
position does not exceed the characteristic size of the 
pocket. In a sense such a retention is analogous to the so 
called cage-effect to be observed for single atoms or small 
molecules in porous glass-like matrices [47-49] . The time 
scale, on which the retention is realized, is limited by 
a time point T esc - This time is a lifetime for the lig- 
and inside the heme pocket in the absence of rcbinding 
and specifies thereby a ligand escape from the pocket. 
Only on a longer time scale (when the ligand succeeds 
to leave the pocket and to migrate over the protein) the 
ligand displacement S(t) starts to increase significantly. 
According to the model, we associate the time r esc with 
the time of particle localization in the hemispheric cage. 
In the simulation the particle displacement S(t) does not 
exceed the cage radius p on the short time scale t < T esc - 
At longer times the quantity S(t) increases with time. 
Due to the diffusion nature of the particle displacement 
the increase in S(t) 2 is a linear function of time and 
S(t) 2 i=a 6Dt at t » T esc . The relation between the time 
T esc and the diffusion coefficient D can be then deter- 
mined from the requirement p 2 ~ S(T esc ) 2 — 6DT esc : 



Tesc — 



6D 



(8) 



Fig. 1 demonstrates a typical temporal dependence 
for the relative particle displacement S(t) 2 / p 2 simulated 
within the framework of our model. The displacement 
S(t) is shown in the figure to be constant (S(t) ~ p) 
at Th <C t <C T esc - On a longer time scale (t 3> r esc ) 
the quantity S(t) approaches asymptotically the diffusion 
law: S(t) 2 w 6Dt — tp 2 /T esc . Notice that for the time 
scale Th <C t <C Tesc the temporal behavior of relative 
displacements S(t)/p is specified by the only parameter 
Tesc- In the following, we will adjust the parameter r esc 
in the simulation. For definiteness, this adjustment will 
be carried out by means of variation in the diffusion co- 
efficient D. The particle mass to, the temperature T and 
the cage radius p will take fixed values typical for the 
ligand and the protein. 

The heme-ligand rcbinding is accepted to be an irre- 
versible process occurring when the ligand is localized 
inside the heme pocket. Therefore this process is sim- 
ulated as a random 'death' for the particle within the 
hemispheric cage. The particle with |r(ifc)| < p is 'oblit- 
erated' if £_k < Ath.k/Treb- Here is a random quantity 
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to be generated for each period Ath,k when the particle 
crosses the cage. The quantity is distributed uniformly 
in the interval [0, 1]. The 'obliterated' particle is excluded 
from the following consideration. 




FIG. 1. A simulated temporal dependence for the 
mean-square relative displacement S(t) 2 /p 2 of particle from 
the initial position at m = 32 amu, T = 300 K, p = 4 A, 
D = 3.2 • 10~ n m 2 /s (solid line). The dashed line shows 
dependence S(t) 2 /p 2 — t/T eac corresponding to the diffusion 
law. The inset shows the quantity Ps/(1 — Ps) as a function 
of the ratio r re b/r esc . (A - our simulation at m = 32 amu, 
T = 300 K, p = 4 A). The solid line gives approximation (9) 
at C s = 0.43. 

The probability P(t) to find the heme in unbound state 
is found as the ensemble-averaged relative number of the 
'non-obliterated' particles at a time instant t. In con- 
trast to the relative displacement S(t)/p, the temporal 
behavior for the probability P(t) depends not only on the 
diffusion properties, but on the rate of hcmc-ligand re- 
binding as well. Hence, the behavior of P(t) can be spec- 
ified in terms of the times r esc and r re b. In general, the 
probability P{t) is a monotonously decreasing function 
of time, which approaches asymptotically a steady value 
P s at t — > oo. This value gives a portion of the hemes 
remaining in unbound state after realization of GR. As 
in diffusion equation (2), in our model the quantity P s 
is a function dependent merely on the ratio between r esc 
and T re ;,. The analysis of simulated data shows that for a 
wide range of values r esc and T re b satisfying the require- 
ment T re b/T esc < 20 (that is, under conditions typical for 
the GR) the best approximation of the dependence can 
be represented by relation 

T^W^cJf^ (9) 

r r s 'esc 

where the coefficient C s is obtained from mean square 
fitting. At to = 32 amu, T = 300AT and p = AA the 
fitting gives a value C s — 0.43. The inset of Fig. 1 
demonstrates a good agreement between approximation 



(9) and the simulated data. 

III. GEMINATE RECOMBINATION OF HUMAN 
HEMOGLOBIN WITH OXYGEN 

We use the described model in order to analyze avail- 
able experimental data for a postphotodissociative reoxy- 
genation of human hemoglobin HbA. The data include 
the measured recombination kinetics and the efficiency 
of oxygen escape from the protein to the solvent for the 
monomolecular (geminate) and bimolecular stages of re- 
combination reaction 

Hb(0 2 ) s + 2 -^Hb(0 2 ) 4 (10) 

going at room temperatures after fast laser-initiated 
breaking of a Fe — 2 bond [51,51]. The kinetic measure- 
ments are carried out with a time resolution of 10 ps for 
a time scale t < 1.5 ns at different fixed pH values for the 
solution. The values of pH fall within an interval between 
6.0 and 9.4: pH = {6.0, 6.8, 7.0, 7.2, 7.7, 8.0, 8.5, 9.4}. 

According to the model, in the analysis of reoxygena- 
tion reaction (10) the temporal decay for the probability 
P(t) is interpreted as a result of two competing processes: 
the heme-oxygen rebinding for the oxygen molecule local- 
ized inside the heme pocket and the diffusion-like migra- 
tion of the oxygen between hemoglobin cavities. Here 
we determine the times T re b and r esc , which specify the 
processes. We determine the times as functions of pH 
and analyze the effect of solution properties on the pro- 
cesses to be considered. Due to the tetramer arrange- 
ment of hemoglobin (the Hb molecule consists of heme 
containing a- and /3-chains) the observed kinetic curve 
represents a reoxygenation kinetics summarized over the 
chains. Here we make no distinction for reaction (10) be- 
tween the a- and /3-chains and determine thereby chain- 
averaged times. 

A. Reoxygenation kinetics for hemoglobin 

The analysis of reoxygenation kinetics for the 
hemoglobin is based on the estimation of the times r re b 
and T esc . The times are found with the help of a nu- 
merical simulation, the iterative procedure for which is 
described above (see Section 2). In the simulation the 
masse of walking particle is accepted to equal the mass 
of oxygen molecule. The temperature T is 300 AT. The 
size p of hemispheric cage is 4 Athat corresponds to the 
time Th < 1 ps. The times T re b and r esc are chosen from 
an interval of values from 0.1 up to 5 ns. The corre- 
lation time r c and the mean length A of free path are 
determined by relations (6) and (7). They are negligibly 
small in comparison with r re b, r esc or p. The simulated 
dependences for the probability P(t) are obtained from 
ensemble averaging for more than 10 6 particles. 
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In the simulation the parameters T re b and r esc are 
so adjusted that the ensemble-averaged temporal depen- 
dence of simulated probability P(t) is the best agreement 
with a measured kinetic curve. The agreement is spec- 
ified by the relative root-mean-square deviation R be- 
tween the simulated and experimental curves. The simu- 
lated dependence for P(t) is shown in Fig. 2 to reproduce 
well kinetic measurements on the considered time scale. 
The minimal deviation R achieved in our calculations for 
each of the fixed pH values does not exceed the mea- 
surement error (R < 0.01). Such an agreement testifies 
that the non-exponential dependence for P(t) with time 
can be explained by a diffusion-like migration of ligand 
over protein matrix. Hence, the parameters T re b and r esc 
can be used for the analysis of the processes, which are 
responsible for the GR. 




f, ns 

FIG. 2. A temporal dependence for the probability P(t) 
to find the heme in unbound state after fast breaking of a 
Fe — O2 bond in human hemoglobin HbA for the geminate 
stage of reaction (10) at pH — 6.0 (solid line - our simulation 
with T re b = 0.495 ns and r esc = 0.84 ns, A - experiment [50]) 
and 8.5 (dashed line - our simulation with T re b = 0.366 ns and 
Tesc = 0.92 ns, o - experiment [50]) . 

The influence of solution properties on the migration 
and the rcbinding of oxygen molecule in hemoglobin is 
assessed from a pH dependence for the obtained times 
T re b and T esc . Our simulation demonstrates a significant 
variation in the rate of hemc-oxygen rebinding with pH 
(see Fig. 3). The increase of quantity pH from 6.0 to 8.5 
results in a shortening for the time T re b by a factor of 1.4 
(from 0.5 down to 0.36 ns). With the following rise of pH 
to 9.4 the parameter T re b appears to increase up to 0.4 ns. 
The minimum magnitude of T re b is observed at pH = 8.5. 
Despite the considerable pH effect for the heme-oxygen 
rebinding, the variation in pH influences slightly on the 
oxygen escape from the heme pocket. The time r esc is 
shown in Fig. 3 to be within a range of values from 
0.82 to 0.92 ns and to be weakly dependent on pH. The 



pii-averaged magnitude of r esc is approximately 0.88 ns. 
Notice that this magnitude is larger than T re b by a factor 
of 2 -=-3. 
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FIG. 3. The times of heme-oxygen rebinding as functions 
oipH for the geminate (r re 6, full circles - our simulation) and 
bimolecular (r s , O _ experiment [50,51]) stages of recombina- 
tion reaction (10). The inset shows the p/f-dependence (A) 
obtained in our simulation for the time r eS c of oxygen escape 
from the heme pocket of HbA. 

The obtained values for the times r re & and T esc are 
in good agreement with the experimental study of the 
alkaline Bohr effect (the variation of the recombination 
rate l/r s with pH for the bimolecular stage of GR) 
[50,51]. The behavior of pii-dependence for the time 
Treb is demonstrated in Fig. 3 to be similar to one for 
the time t s of bimolecular rebinding. Such an agree- 
ment testifies that for the monomolccular GR the vari- 
ation of the rcbinding rate with pH can be associated 
with the same structural transformation as for the bi- 
molecular stage of reaction (10). Histidine imidazoles of 
C-terminal sites and a- amides of TV-terminal sites seem 
to be the aminoacid residues, which are responsible for 
this transformation [2,8,53]. Specifically, in the alkaline 
Bohr effect the interaction between the solvent and the 
pi46His residue (a C-terminal histidine of /3-chain) is 
one of the most probable reasons for the heme struc- 
ture modification and the rearrangement of neighboring 
aminoacid residues [2] . Our simulation confirms that the 
variation in pH can result in essential structural trans- 
formations in immediate proximity from the iron atom. 
The strong pii-dependence for the times T re b and t s is a 
consequence of the transformations. 

The ligand penetration from the solvent into the heme 
pocket is shown for 2 or NO in Hb [5] to be a process 
restraining the rate of bimolecular recombination (10). 
Therefore, the similarity between the piJ-dependences 
for Treb and t s testifies that the change in pH has a 
slight effect on the oxygen migration in hemoglobin at 
the mono- and bimolecular stages of recombination (10). 
The weak pii-dependence for the obtained times of oxy- 
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gen escape from the pocket confirms this assumption. 
Such a pfl-invariant behavior for the oxygen migration 
can be interpreted by the independence of mobility for 
the hemoglobin side chains (which seem to be responsi- 
ble for the ligand transitions between cavities of heme 
protein [38]) on pH of the solution. 



B. Efficiency of oxygen escaping from hemoglobin 

The obtained times r re b and T esc are used then in 
order to estimate the efficiency of oxygen escape from 
hemoglobin as a function of pH . The efficiency is pro- 
portional to the quantum yield of photodissociation and 
can be associated with the portion P s of the hemes re- 
maining in unbound state after realization of the gemi- 
nate reoxygenation stage [54]. In our simulation the ra- 
tio T re b/T esc falls within a range of values from 0.4 up 
to 0.6. It implies that the portion P s can be determined 
from approximation(9). The behavior of p_ff-dependcncc 
for the obtained quantity P s agrees well with one for 
the measured quantum yield of photodissociation [51,51]. 
The quantity P s is shown in Fig. 4 to be proportional to 
the apparent quantum yield 7 for the whole investigated 
scale of pH . 
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FIG. 4. The portion P a of the hemes remaining in unbound 
state after realization of the geminate stage of HbA reoxy- 
genation (10) (full circles - our simulation) and the apparent 
quantum yield 7 of HbA photodissociation (O - experiment 
[50,51]) depending on pH values of the solution. 

Notice that for the studied range of pH values the 
quantity C s T re b/T esc is considerably low in comparison 
with 1 {C s T re b I r esc ~ 0.15 -j- 0.25) and the time T esc is 
practically constant. Therefore the relation of the por- 
tion P s with the time r re b is close to the linear law: 



s 7"re&- 
7~esc 



(11) 



In our study (see Fig. 3 and Fig. 4) the p-ff-dependences 
obtained for P s and similar that testifies again 



that the transport properties for the oxygen molecule in 
hemoglobin do not depend on pH . 



C. Diffusion properties of oxygen migration in 
hemoglobin 

The analysis of X-ray diffraction data [55,56] for oxy- 
genated and deoxygenated species of human hemoglobin 
(PDB ID 1HHO and 2HHB, correspondingly) shows 
that the cage radius p to be associated with the heme 
pocket size is a quantity ranging from 1 up to 5 A (taking 
into account the Van der Waals radiuses). Hence, 
the diffusion coefficient D for the oxygen migration in 
hemoglobin can be estimated from relation (8): D = 



p 2 /6i 



0.2 -f- 5 • 10 11 m 2 /s. This coefficient is 



termccliatc to diffusion coefficients for small molecules in 
water ( 10~ 9 m 2 /s) and solids (10~ 18 m 2 /s at T < A00K) 
[57] 

According to the diffusion law, at a time instant t m the 

root-mean-square displacement \J (|r(t)| 2 ) of the ligand 

from the iron is approximately equal to p\/t m /T esc . It 
implies that on the completion of kinetic measurements 
(t rn = 1.5 ns [50]) the oxygen remains inside the protein 
and is localized in immediate proximity from the heme 

pocket: \J (|r(t)| 2 ) 1.3p < 7 A. This conclusion is 
consistent with results of spectroscopy investigation of 
motional dynamics for CO in Hb [42]. 



IV. CONCLUSION 

We have represented a simple model of the geminate 
recombination of heme proteins with small ligands. The 
model takes into account dynamic properties of ligand 
displacement in protein matrix. In the model the recom- 
bination is due both to the heme-ligand rebinding and 
to the diffusion-like migration of ligand between protein 
cavities. The temporal behavior for the probability P(t) 
to find the heme in unbound state is specified in terms of 
two parameters. They are the time r re b of heme-ligand 
rebinding for the ligand inside the heme pocket and the 
time T esc of ligand escape from the pocket. 

We have applied our model in order to analyze a 
postphotodissociative geminate reoxygenation of human 
hemoglobin at various pH values of the solution. The 
measured kinetic curves and the efficiency of oxygen es- 
cape from the hemoglobin are well reproduced in our sim- 
ulation. It testifies that the non-exponential behavior for 
the probability P(t) can be explained by a diffusion-like 
migration of ligand over protein cavities. This conclu- 
sion is consistent with recent kinetic measurements [58]. 
We believe that the theory-experiment agreement may be 
considered as an additional validation for the glass-like 
model of proteins. 
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Our study demonstrates also that the variation in pH 
can result in considerable changes for the rate of hcme- 
ligand rebinding. At the time, the oxygen migration in 
hemoglobin depends slightly on pH . We have interpreted 
this effect as a result of essential structural transforma- 
tions in immediate proximity from the iron atom. Cer- 
tainly this conclusion demands a more detailed and thor- 
ough examination. In any case we suppose that the pH- 
induced modification of the initial stage of GR (if the 
modifications are observed) can be explained by a change 
in the rate of heme-ligand rebinding. 
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